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Abstract: In this paper we show that there are applications that transform the 
movement of a pendulum into movements in R^. This can be done using Euler top 
system of differential equations. On the constant level surfaces, Euler top system 
reduces to the equation of a pendulum. Those properties are also considered in the 
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case. Another aspect presented here is stochastic Euler top system of differential 
equations and stochastic pendulum. 
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1 Introduction 

The dynamics of some mechanical systems is described using the rigid body dy- 
namics with a fixed point, mathematical pendulum or oscillators. These systems 
belong to a class of differential equations from R.^ with the right side polynomial 
functions of degree greater or equal to two. From this category we will consider 
Euler top system of differential equations. We begin our study from mathematical 
pendulum (and its variants: with delay, fractional and stochastic) approach. 

The Euler top system of fractional differential equation belongs to a class of 
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differential equations that are described using polynomial functions. It has the form 

' Xi{t) = X2(t)X3it), 

• m) = -xi(t)x3(t), (1) 

Mt) = Xi(t)X2it). 

Because system ^ has three Hamilton-Poisson realizations, three conservation 
laws are given by the Hamiltonians Hi, H2 and [i2|: 

1. Hiixiit),X2it),X3it)) := + xjit)); 

2. H2(xi(t), X2(t), X3(t)) := -^(xjit) + xjit)); 

3. H^iXiit), X2(t), X3{t)) ■= x^it) - xjit), 

and the other three conservation laws are given by the corresponding Casimir func- 
tions of the above realizations 

1. Ci(xi(t),X2(t),X3(t)) := ^(xjit) + xjit)); 

2. C2(xi(t),X2(t),X3(t)) := ^(x^^(t) - xjit)); 

3. C3(xi(0,X2(0,-^3(0) := xl(t) + xjit). 

A simple mathematical pendulum is the mathematical model of a ball, having 
the mass m, which hangs in a point O by a bar of length /, and the point O performs 
movement in a plane [9|. 

The Euler-Lagrange equation that describes the movement of a pendulum is 
given by 

I6(t) + g sin Bit) + xq cos 6{t) - yo sin e(t) = 0. (2) 
The dumping pendulum equations with periodic force is 

N 

m + 2h sin eit) + cos e{t) + f2{t) sin e{t) + ^ ape{t)\e{t)r^ = 0, (3) 

and for /i := 0, /2 := and ap := 0, p = I...N, then ® reduces to 0{t)+2h sin 0(0 = 
0. 

In the first section we will determine the analytical solutions for Euler top sys- 
tem taking into consideration the conservation laws that it owns, and point out the 
analytical solution for pendulum. In the second section we have presented the co- 
nnection between Euler top system and pendulum: the restriction of the system to a 
constant level surface represents the pendulum equations. The third section presents 
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the Euler top system of differential equations with delay argument, along the OZ and 
OX axes. These new systems have also conservation laws and the restriction of the 
orbits at these surfaces of constant level determined by the conservation laws are 
mathematical pendulums with delay argument. In the forth section we presented 
the Euler top system of fractional differential equations. We have used Caputo frac- 
tional derivative in OZ and OX directions. As in the previous case, this system of 
fractional differential equations have conservation laws and the restriction of the 
system to the constant level surfaces is a fractional pendulum. In Section 5 we pre- 
sented stochastic Euler top system and stochastic pendulum. We considered Ito and 
Stratonovich integrals for describing the stochastic process, using a Wiener pro- 
cess. For all these cases numerical simulations are done. In the last section some 
conclusions are presented and ideas for future work. 

2 Euler top system and simple pendulum - analytical 
solutions 

Let us consider the Euler top system of differential equations O and the integrals 
of motion given by 

xlit) + xjit) = 2H^, xlit) + xlit) = 2K^. (4) 
From @, results that 

x](t) = 2H^ - xl(t), xlit) = 2K^ - xlit). (5) 
Replacing ^ in the first equation in ([T]) we get: 

{X2)\t) = (x,f + (xsfit) = {2H^ - xlm2K^ - xi(t)) (6) 

and so, 

t = I — — du. (7) 

that shows that X2{t), respectively xi{t) and X2{t) are elliptic functions of time [|8]|. 

In the case when H = K, the quartic under the square root has double roots and 
dV]) can be explicitly integrated by means of elementary functions in the following 
manner. The equation 

X2{t) = +{2H^ - xim 
with X2(0) = 0, has the solution 

X2{t) = +H V2 tanh(/f ^f2t). (8) 
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Substituting ^ in ©, we get 

jci (t) = +// V2 sech(// V20, X3(0 = +// V2 sech(// a^O- (9) 

So, the equations dH) and ([3]) represent the two heteroclinic orbits for the Euler top 
system and are given by 

( + // V2 secUH V20, +H V2 tanh(// V20, +// V2 sech(// V20). 

In the case when H K, the integral ^ can be computed using Jacobi's elliptic 
functions [|7]|. We use relations 

— sn M = cn M dn M, cn u = \ - sn u, dn w = 1 - m sn m 
at 

and 

X2(t) = // V2 sn (// ^^2t■ -^), (10) 

with the initial condition X2(0) = 0. Choosing the time deviation, appropriately, we 
can assume that ^2(0) > 0- From Q results that 

a/77 a/77 
jci(0 = //V2cn(//V2f;-— ), = K^^niHyflt;-—]. (11) 

If (f) denotes the period invariant of Jacobi's elliptic functions, then x\{t) and X2{t) 
have the period A^jH y/l, whereas ^3(0 has the period 20/// V2. 

Proposition 1 a) If H = K, than Euler top system ^ has an analytical solution 
given by © and ®; 

b) IfH K, than the Euler top system has the analytical solution given by (flOl) 
and (fTTT) . 

□ 

Proposition 2 [1] The analytical solution for simple pendulum 9{t)+2h sin 6(t) = 0, 
with initial conditions 9(0) = 6q and 9(0) = is given by 

9(t) = 2 arcsin j sin 9o sn ( sin^ y ) - ^/2ht■, sin^ y }. 
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3 Euler top system and simple pendulum 

In this section we will show the way the Euler top system and the simple pendulum 
are linked. We will show that the movement of the Euler top system is reduced 
to pendulum movement on the constant level surfaces H and K, described by the 
conservation laws: 

\ixi{t)f + \{x2{t)f = H, (12) 

\(X2(t)f + ^(X3(t)f = K. (13) 

Since H and K are conserved, the Euler top motion takes place along the inter- 
sections of the level surfaces of the energy and the angular momentum in M?. 

Proposition 3 Let us consider the Euler top system of differential equations O. 

1. The function H, given by (fT2l) is a conservation law for system ([T]); 

2. The solution of ([U) on the constant level surface defined by (fT2]) . given by 

(xiit)f + (X2(t))^ = 2H = const, H>0 (14) 

is 

x,{t) = V2HCOS X2it) = V2Hsin^, ^3(0 = -^m (15) 

where 6(t) is the solution of pendulum equation 6{t) + 2H sin 0(f) = 0. 
Proof: 

1. By deriving (fT^ and by replacing it in (O, we have 

H(t) = Xiit)Xi(t) + X2it)X2it) = 0. 

And so is a conservation law. 

2. Using a direct calculus, it can easily checked that (fT5l) is a solution for ([T]) 
and reciprocal. □ 

Proposition 4 1. The function K, given by (fT3l) z^' a conservation law for the 
Euler top system (O; 
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2. The solution of on the constant level surface defined by (fT3l) . given by 

(X2(t)f + {X3(t)f = 2K = const, K>0 (16) 

is 

1 . I — eit) I — e(t) 

Mt) = -^e{t\ X2{t) = V2ZC0S ^3(0 = V2Zsin (17) 

where 6{t) is the solution of pendulum equation 6(t) + 2Ksm6(t) = 0. 
Proof: 

1. By deriving (fT3l) and by replacing it in ([U), we have that ^ is a conservation law 
because 

K(t) = X2{t)X2(t) + x^{t)x^{t) = 0. 

2. By direct calculations, it can be easily checked that (fTTl) is a solution for ([T]) and 
reciprocal. □ 

Remark 5 The dynamics of Euler top system of differential equations in is a 
union of two-dimensional simple pendula. 

□ 

For the initial conditions Xi(0) = 0.1, X2{0) = 0.1 and .^3(0) = 0.2, the Euler 
top system is represented in the first figure and the pendulum is represented for the 
initial condition 0(0) = -3.8. 



EULER TOP SYSTEMOFFRACTIONALDIFFERENTIAL EQUATIONS PENDULUM REPRESENTATION 
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4 Euler top system and simple pendulum - with delay 
argument and fractional derivative 

A differential equation with delay argument is defined in [|4||. A second order diffe- 
rential equation with delay argument is given by 

m = csm(eit-T)), (18) 

where c € R is a solution of a differential equation on the circle 5 ^ = {y 6 R^|jj + 
yl = \},6 an angle variable determined up to a multiple of 2n, and t > 0. 

From Proposition [3] and Proposition!?] we can deduce the following results. 

Proposition 6 The Euler top system of differential equations with delay argument 
is given by 

' X-xit) = X2(.t)X3(t), 

< Mt) = -Xiit)X3(t), (19) 
^ X^it) = Xi{t-T)X2(t-T). 

The system (fT^l ) has the following properties 

a) The function H given by (fT2)) ; 

b) The solution of system ([19] ) on the constant level surface ([T2]) is given by (fT5]) . 

where 6{t) is the solution of 

e{t) + 2Hsme(t-T) = (20) 

and reciprocal. □ 

Proposition 7 The Euler top system of differential equations with delay argument 
given by 

' Xyif) = X2{t - T)X3(t - t), 
• X2(t) = -Xi{t)X3(t), (21) 
. i3(0 = Xiit)X2(t), 

has the following properties 

a) The function K given by (fT3]) is a conservation law for the system (|2T]) .- 

b) The solution of system (|2T]) on the constant level surface ([T3]) is given by ([T7]) . 

where 6{t) is solution of the equation 

e(t) + 2Ksm9{t-T) = 0. (22) 



8 



□ 

This system is considered to be a starting point in studying differential equations 
with delay argument for differential manifold. 

For H = 0.5, and r = 1, the pendulum equation with delay argument and with 
initial condition 6(0) = 2, is represented in the following figure. The Euler top 
system with delay argument (fT9l ) is represented in the second figure, for the initial 
conditions Xi(0) = 0.1, ^2(0) = 0.05, ^3(0) = 0.2. 




.4 0.5 

X3 



Pendulum with delay r = 1, // = 0.5 Euler top system with delay r = I H = 0.5 

For K = 0.3, and r = 1, the pendulum equation with delay argument and with 
initial condition 6(0) = 2, is represented in the following figure. The Euler top 
system with delay argument ([T3 is represented in the second figure, for the initial 
conditions ;ci(0) = 0.1, X2(0) = 0.05, ^3(0) = 0.2. 




Pendulum with delay t = 1, ^ = 0.3 Euler top system with delay t = I K = 0.3 
Using Caputo fractional derivative fT\, the following propositions take place. 
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Proposition 8 The Euler top system of fractional differential equations, given by 

' y-xif) = X2(t)X3(t), 

X2(t) = -Xi(t)X3(t), (23) 
[ D^X^it) = Xiit)X2it), 

with a 6 (0, 1), has the following properties 

a) The function H is a conservation law for (|23]) ; 



b) The solution of the system (|23l) on the constant level surface (fT2)) . with 6{t) is the 
solution of the fractional equation 

D'^^^e(t) + 2H sin e(t) = 0, (24) 

and reciprocal. 

□ 

Proposition 9 The Euler top system of fractional differential equations, given by 

D"x,{t) = X2(t)X3(t), 

Ht) = -xi{t)x3(t), (25) 

X3(t) = Xi(t)X2(t), 

with a 6 (0, 1), has the following properties 

a) The function H is a conservation law for (|25l) ; 



b) The solution of the system (|25]) on the constant level surface ([13]) . with 0(t) is the 
solution of the fractional equation 

D'^^^eit) + IK sin e{t) = 0, (26) 

and reciprocal. 

□ 

By using the Adams-Moulton method for integration, for the initial condition 
9{0) = -3.1, the solution of the fractional differential equation (|24l) is represented 
in the following graphics for a = 0.8, respectively for or = 1. 
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FHiCTIONAL PENDULUM(alpha 0.8, 



FRACTIONAL PENDULUM (alpha- 1) 

r\. n. n, n 
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It can be observed that the pendulum solution is asymptotically stable for < 
a < I and it is oscillatory for or = 1. 

The solution for the system of fractional differential equations (1231) . respectively 
for (1251) . is represented in the above graphics, for the initial conditions Xi(0) = 
0.1, X2i0) = 0.1 and XsiO) = 0.3. The cases of a = 0.8 and a = I are illustrated. 




EULER TOP SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS (alpha=0.8) EULER TOP SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS (alpha=1) 
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5 Stochastic Euler top system and stochastic pendu- 
lum 

A Wiener process describes rapidly fluctuating random phenomena. Stochastic di- 
fferential equations (SDE) are stochastic integral equations and are written symbo- 
lically in a difl'erential form. We will consider such a Wiener process of the form 

dx{t) = f{x(t))dt + g{x{t))dW{t), (27) 

where / is the slowly varying continuous component called drift coefficient and g is 
the rapidly varying continuous component called diffusion coefficient. The integral 
representation is of the form 

x{t) = x{to)+ r f{x{s))ds+ r g{x{s))dW{sX (28) 

where W{t)h a Wiener process, a Gaussian process with W{0) = and A^(0, ^-distributed 
W{t) for each ? > 0, so 

E(iy(0) = 0, ^{{W{t)f) = t. 

The first integral is a Riemann-Stieltjes integral and the second one is a stochas- 
tic integral. The most studied interpretation of the stochastic integral are those of 
Ito and Stratonovich. The choice of interpretation depends on the type of analysis 
required for solution dH. Ito stochastic calculus is closely related to diffusion pro- 
cesses and martingale theory [0. The solution of (ITTl) is a diffusion process with 
transition probability p = u{x{t)), satisfying the Fokker-Planck equation 

d d 1 

^u{x{t)) = --—f{x{t))u{x{t)) + -— — [(g(x(0)/(x(0)w(40)]. (29) 
ot ox{t) 2 d{x{t)y 

Equations (1271) and (l29l) contain the same statistical information from a one-particle 
process point of view (but not if we think the Ito equation as describing a random 
dynamical system) Q. 

An Ito SDE is written in the form (l27l) and a Stratonovich SDE is written sym- 
bolically in the form 

dx{t) = f(x{t))dt + g{x(t)) o dW(t), (30) 
and in the integral form as 

x(t) = x(to)+ f f(x(s))ds+ f g(x(s))odW(s). (31) 

Jto Jto 
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It is possible to switch between these two approaches, in the sense that the Ito 
SDE (ETl) has the same solution as Stratonovich SDE 

dx(t) = f{x(t))dt + g{x{t)) o dW{t), (32) 

with modified drift coefficient 

mm = mm - l.g(x(m^Ax(m- 

— 2 dx{t) 

If W{1), W(d) are d independent Wiener processes, and x{t) = {xi{t), Jc„(0) 
then the multi-Wiener process case can be written in the form 



dxi(t) = fiximdt + J]gij(x{mdwu), 



(33) 



7=1 



with g{x{m anx d matrix and dW adx \ matrix. 

In Stratonovich case, the stochastic system of differential equations with a multi- 
Wiener process, can be written in the following manner 



dxit) = mmdt + gjixim ° dwu). 



(34) 



where 



^ n d 



mm = f(x(m -2LL ^'^^^^^^^^'^ 



it) 



k=i j=i 

The Euler top system of stochastic differential equations can be represented in 
the following form. 



' dxi(t) = X2{t)x^{t)dt + xi{t)dW^{t) 
dx2{t) = -xi{t)x^{t)dt, 
dx-iit) = Xi{t)x2(t)dt + dW^{t), 



(35) 



with the Wiener process W(t) = (W\t),0, W^(m, the drift coefficients f\x{t)) 

X2(t)x3it),fix{m = -xiit)x3{m f\x{m = ^1(0^2(0, ^(o = (xi(o,^2(o,^3(0)^, 

f{x{m = {f^{x{m^P{x{t)),f{x{mf and the diffusion coefficient vectors 





( x,{t) \ 




{ ^ 1 




{ ^ \ 


\x{m = 





, g\x{m = 





, g\x{m = 







I J 




lo J 




1 1 J 
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The corresponding (Ito) Fokker-Planck equation for the probability density p = 
u{x{t)) reads 

d d d d 

— U(x(t)) = -^—[X2(t)X3(t)u(xit))] + ——[Xy(t)X3it)uix(t))] - ^—[Xyit)X2(t)u(x(t))] 
Ot OXi{t) OX2(t) OXjit) 



1 



+ -- 



2d(xm- 



■[{xmMxm + ;7 



1 



2dix,it)) 



:u{x{t)). 



In the Stratonovich case, stochastic system (1351) can be written using relation 
(|34l) in the following manner 



dxiit) = {X2{t)x^{t) - \x2{t))dt + Xi{t)dW^{t) 
dx2(t) = -Xi{t)x3{t)dt, 
dx^it) = Xi(t)x2(t)dt + dW^it). 



(36) 



The stochastic system (|35I) . respectively (|36l) . is implemented in Matlab, using 
Milstein scheme, for initial conditions xi{l) = 0.1, ^2(1) = 0.1, ^3(1) = 0.1 and 
orbits are represented in the following figures. 
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0.1 0.3 0.3 
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If the SDE of Euler top system has the form 



' dxiit) = X2(t)x3(t)dt + y/xi(t)dW\t) 
dx2{t) = -xi(t)x^(t)dt + y/x2(t)dW\t), 
dxsit) = Xi(t)x2(t)dt + y/x3(t)dW^it), 



(37) 



then drift coefficients are 



fiit)) = X2(t)X2it), f(xit)) = -Xi(t)X3it), f(x(t)) = Xi(t)X2it), 

withxW = (^,(0,^2(0,^3(0)^ mt)) = {f{x{t))J\x{t))J\xm\ and the dif- 
fusion coefficient vectors 





( y/xiit) ] 




f ^ ] 




f ^ \ 


\x{t)) = 





, ^'(-^(0) = 


ylx2{t) 


, ^'(-^(0) = 







I J 




[ J 




y ^Ix^it) j 



then the associated (Ito) Fokker-Planck equation for the probability density p 
u(x(t)) is 



d 

—u{x{t)) 
ot 



d d d 

[x2(0-^3(0w(^(0)] + Tr^[^i(0^3(0w(-«(0)] - Tr^[-^i(0-^2(0"(-^(0)] 



5^1 (0 5^2(0 5^3(0 

1 52 1 52 1 52 

+ . ..A xi{t)u{x{m + ,.,,A x2{t)u{x{m + - 



2d{x,{t)f 25(X2(0)' 25(^3(0)" 

The Stratonovich stochastic Euler top system is written in the following way 



[x3(0m(^(0)]- 



dxi{t) = {X2{t)x2it) - \)dt + y/xi{t)dW\t) 
dx2(t) = -(xi(t)x3it) + \)dt+ y/x2it)dW^{t), 
dx^it) = (xi(t)x2(t)dt + |) + y/xi(t)dW^(t), 



(38) 



Stochastic system dJT]) . respectively (l38l) . can be implemented using stochastic 
Euler method which represents a square-root model. For initial values Xi{l) = 
1, JC2(1) = 0.8, ^3(1) = 0.2, orbits are represented in the following figures. 



y 



m 
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The stochastic pendulum equation is considered in the following manner. The 
dynamics of a non-dissipative classical pendulum of the form 6{t) + 2H sin d{t) = 0, 
can be expressed as a system of stochastic differential equations expressed like 

dxi(t) = X2(t)dt+ y/xi(t)dW\t), 
dx2(t) = -2H sm(xi(t))dt + yfx^dW^(t), 

and the Stratonovich stochastic pendulum equations are 

dxi(t) = (X2(t) - \)dt + y/xi(t)dW\t), 
dx2(t) = -(2H sm{xi(t)) + \)dt + ^[x^dW'^{t), 

For the probability density p = u{x{t)), the corresponding (Ito) Fokker-Planck 
equation is given by 



(39) 



(40) 



d 



d d 

{X2{t)u{x{ty)\ + -——{2H?,m{xi{t))u{x{t))] 



+ o ^r^;^[^i^^)"WO)] + ,^..2 ^x2it)u(xit))]. 

2d(xi(t)y 2d{x2{t)y 

Using stochastic Euler method on square root process, for initial conditions 
x\{\) = 1, ^2(1) = 0.8 we get the following graphics for stochastic systems (1^ 
and (go]) 



0.1 0.2 0.3 
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6 Conclusions 

In this paper we presented the Euler top system in and the mathematical pen- 
dulum, but also the connections between them: the existence of some applications 
that transform the movement of a pendulum into a movement in M?. That means that 
the restriction of the Euler top system on a constant level surface is the pendulum 
equation. This property is also true in the case of Euler top system of differential 
equations with delay argument, respectively mathematical pendulum with delay ar- 
gument, and in the case of fractional system of differential equations, respectively 
fractional pendulum. We have also studied the Euler top system and mathematical 
pendulum from the stochastic point of view, using Ito and Stratonovich integrals 
for a Wiener process. Numerical simulations were done using Maple 12 and Mat- 
lab. In the case of fractional Euler top system and fractional pendulum we used the 
Adams-Moulton integration method for their representation, and in the stochastic 
case we used the Milstein scheme, that is a convergent numerical algorithm. In the 
future we will study other aspects of these problems, such as stochastic Lyapunov 
functions, stochastic Lyapunov exponents for determining the stochastic stability in 
the equilibrium points of a considered system, classical, with delay of fractional. 
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